DETERMINING AN UNBOUNDED POTENTIAL 
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GEOMETRIES 
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Abstract. In anisotropic inverse problems were considered 
in certain admissible geometries, that is, on compact Riemannian 
manifolds with boundary which are conformally embedded in a 
product of the Euclidean line and a simple manifold. In particular, 
it was proved that a bounded smooth potential in a Schrodinger 
equation was uniquely determined by the Dirichlet-to-Neumann 
map in dimensions n > 3. In this article we extend this result to 
the case of unbounded potentials, namely those in L"/^. In the 
process, we derive Carleman estimates with limiting Carleman 
weights similar to the Euclidean estimates of Jerison-Kenig [lO] 
and Kenig-Ruiz-Sogge [H]. 
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1. Introduction 

In this paper we consider the problem of proving U' estimates for 
limiting Carleman weights on Riemannian manifolds, and the related 
inverse problem of recovering an L"/^ potential from the Dirichlet-to- 
Neumann map (DN map) related to the Schrodinger equation. The 
main motivation comes from the inverse conductivity problem posed by 
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Calderon [T]. This problem asks to determine the conductivity function 
of a medium from electrical measurements made on its boundary. 

In mathematical terms, if d R" is the medium of interest hav- 
ing a positive conductivity coefficient 7, in the Calderon problem one 
considers the conductivity equation 

V ■ 7VM = inn 

and defines the DN map by 

du 

: u\dn h-^ 7— 

ou an 

This operator maps the voltage at the boundary to the current given 
by 7 times the normal derivative, which encodes the electrical mea- 
surements at the boundary. The inverse problem of Calderon asks to 
determine 7 from the knowledge of A^. This problem has been exten- 
sively studied and we refer to [50] for a recent survey. 

The anisotropic Calderon problem considers the case where the con- 
ductivity 7 is a symmetric positive definite matrix instead of a scalar 
function. This corresponds to situations where the electrical properties 
of the medium depend on direction. The problem is open in general in 
dimensions n > 3, see |1] for known results and more details. Follow- 
ing [18] the problem may be recast as the determination of the metric 
g on a compact Riemannian manifold (M, g) with boundary from the 
corresponding DN map. In [4j progress was made on the anisotropic 
Calderon problem in the following class of conformal smooth manifolds. 

Definition. A compact Riemannian manifold {M,g), with dimension 
n > 3 and with boundary dM , is called admissible z/ M d R x Mq for 
some {n — 1)- dimensional simple manifold {MQ,go), and if g = c{eQ)go) 
where e is the Euclidean metric on R and c is a smooth positive function 
on M. 

Here, a compact manifold {Mq, g^) with boundary is simple if for any 
p e Mq the exponential map exp^ with its maximal domain of definition 
is a diffeomorphism onto Mq, and if OMq is strictly convex (that is, the 
second fundamental form of OMq ^ Mq is positive definite). 

In [1] it was proved that a Riemannian metric in a conformal class 
of admissible geometries is uniquely determined by the DN map. This 
was obtained as a corollary of a result for the Schrodinger equation in 
a fixed admissible manifold, stating that a bounded smooth potential q 
is determined by the corresponding DN map. In [3] all coefficients were 
assumed infinitely differentiable. In this paper we relax this require- 
ment and show that a complex potential q e L"'/^(M) is determined 
by the DN map. 
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To state the main result, assume that (M, g) is a compact Riemann- 
ian manifold with smooth boundary dM, and let Ag be the Laplace- 
Beltrami operator. Given a complex function q E L"/^(Af), where 
n > 3 is the dimension of the manifold M, we consider the Dirichlet 
problem 

(-Ag + q)u = in M, uIqm = /• 

We assume throughout that is not a Dirichlet eigenvalue for this 
problem, and then standard arguments (see Appendix |Al) show that 
there is a well-defined DN map 

A,,, : H'/\dM) ^ H-'/^idM), f ^ O^uUm- 

The following uniqueness theorem is the main result for the inverse 
problem. (The assumption q G L"/^ may be considered optimal in the 
context of the standard wellposedness theory for the Dirichlet problem 
with potentials, and also for the strong unique continuation principle 
to hold [lOj.) 

Theorem 1.1. Let {M,g) be admissible and let qi,q2 be complex func- 
tions in U'-/'^{M). If Ag^q-^ = Ag g^, then qi = q2- 

In the case where M is a bounded domain in R"' and g is the Eu- 
clidean metric, this result is due to Lavine and Nachman [TT] following 
the earlier result of Jerison and Kenig for qj G L"/^"'"'^(M) for some 
e > (see Chanillo [2] for an account and also for a related result with 
qj in a Fefferman- Phong class with small norm). As mentioned above, 
if g is a smooth function on an admissible manifold M this result was 
proved in [3] by using Carleman estimates. In fact, smoothness of q 
is not essential, and by inspecting the proof of ^ the uniqueness result 
can be extended to bounded continuous q (with the complex geomet- 
rical optics construction in the proof going through for q G L"-[M)). 
However, the proof for q G L*^/^ requires to replace the Carleman 
estimates in |3] with corresponding Carleman estimates. 

The other main result in this paper is a L'^ Carleman estimate for 
limiting Carleman weights on Riemannian manifolds. The concept of 
limiting Carleman weights was introduced in [T^ as part of a general 
procedure for producing special complex geometrical optics solutions to 
elliptic equations, with applications to inverse problems. We refer to [3] 
for a precise definition and more careful analysis of limiting Carleman 
weights, also on Riemannian manifolds. For present purposes it is 
sufficient to mention that the existence of a limiting Carleman weight 
on (M, g) in dimensions n > 3 is locally equivalent with the manifold 
being admissible, and that typical limiting Carleman weights in R", 
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n > 3, include the linear weight <f{x) = xi and logarithmic weight 

(p{x) = log 

The last two weights are featured in the literature of Carleman esti- 
mates and unique continuation, in particular in the scale invariant 
Carleman estimates of Kenig-Ruiz-Sogge ^12] for the linear weight and 
of Jerison-Kenig [10] for the logarithmic weight. We prove an analogue 
of these estimates on more general Riemannian manifolds. Note that 
the existence of a limiting Carleman weight requires at least locally a 
product structure on the manifold, and therefore the following result is 
stated for the linear weight on a product manifold. The result, in the 
case when the manifold (Mq, go) below is the standard n—1 dimensional 
torus, is due to Shen [231. 



Theorem 1.2. Let {MQ,go) be an {n — 1)- dimensional compact man- 
ifold without boundary, and equip R x Mq with the metric g = e (B Qo 
where e is the Euclidean metric. The Euclidean coordinate is denoted 
by Xi . For any compact interval ICR there exists a constant C/ > 
such that if \t\ > 4 and 



t' i Spec(-A 



then we have 



In 



'L7^{RxMo) — " ^ "L7JT5(RxMo) 

when u G C^[l x Mq). 

The proof of the Carleman estimates for limiting Carleman weights 
in [3] is based on integration by parts and cannot be used in the U' 
setting. However, in [I3j another proof of the Carleman estimate is 
given; this proof is based on Fourier analysis and gives an explicit in- 
verse for the conjugated Laplacian. We will derive the LP bounds from 
this explicit inverse operator. This follows the proof of the Carle- 
man estimate of Jerison-Kenig [10] using Jerison's approach [9], [261 
Section 5.1] based on the spectral cluster estimates of Sogge [26]. Fi- 
nally, if one allows strongly pseudoconvex Carleman weights then much 
stronger estimates are available (see for instance [l5l[T6]), however for 
the applications to inverse problems it seems necessary to restrict to 
limiting Carleman weights. 

Remark 1.3. The above theorems are in the setting of (conformal) 
product manifolds. However, the results also apply to warped products. 
If / : R — )■ R is a smooth function and (Mq, go) is an (n— l)-dimensional 
manifold, the warped product R Xg2/ Mq is the manifold M = R x Mq 
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endowed with the metric 

We choose coordinates yi = ri{xi), y' = x' for a suitable smooth strictly 
increasing function r/. In fact, if 

r^(t) = /" e-^(^) ds 

then r]'{t)~'^ = e^-^*^*^ and the metric in y coordinates becomes a confor- 
mal multiple of a product metric, 

Warped products have a natural limiting Carleman weight (p{y) = yi, 
and Theorem 11.11 remains true in conformal multiples of warped prod- 
ucts whenever (Mo,(?o) is a simple manifold. 

The paper is organized as follows. Section [1] is the introduction. In 
Section [2] we prove the Carleman estimate complemented with the 
usual Carleman estimates. Section |3] presents the construction of 
complex geometrical optics solutions for Schrodinger equations with 
L"'^^ potentials in admissible geometries. The proof of Theorem 11.11 is 
contained in Section HJ modulo a uniqueness result for an analogue of 
the attenuated geodesic ray transform acting on unbounded functions. 
This last result has a different character than the rest of the proof, 
and it is therefore established separately in Section |5l There are two 
appendices concerning the wellposedness of the Dirichlet problem and 
the normal operator for the attenuated ray transform. 

Acknowledgements. The last named author would like to thank 
Adrian Nachman for explaining his unpublished argument with Richard 
Lavine [17] which proves a uniqueness result for L"/^ potentials in Eu- 
clidean space. C.K. is supported partly by NSF grant DMS-0968472, 
and M.S. is supported in part by the Academy of Finland. D. DSF. 
would like to thank the University of Chicago for its hospitality. 

2. LP Carleman estimates 

The aim of this section is to prove Theorem II. 2[ which is an analogue 
of the LP Carleman estimates obtained in the Euclidean case by Jerison 
and Kenig [10] (for logarithmic weights) or by Kenig, Ruiz and Sogge 
[T2] (for linear weights). In fact, we prove a more general result which 
implies Theorem 11.21 by taking / = e™^ A^e'^^^^w for u G C^{I x Mq). 



The case when (Mcf^o) is the standard n — 1 dimensional torus is due 
to Shen 



Proposition 2.1. Let I (1 H be a compact interval and (Mo, (70) a 
compact {n — 1)- dimensional manifold without boundary. Equip N = 
I X Mo with the product metric g = e Q) go- For |r| > 4 with ^ 
Spec{—Ag^) , there is a linear operator Gr '■ L'^{N) — )■ H'^{N) such that 

e™i(-Ag)e-^^-iG^t; = v for v E L^{N), 

G,e™i(-A3)e-™-it; = v for v e C^{N''''). 

This operator satisfies 

Co 

WGrfWh^iN) < |;^||/||l2(7V), 
l|G'r/||//i(Ar) < C'o||/||l2(^-), 

^ < Coll/ll ^ , 
where Cq is independent of t (but may depend on I). 

Remark 2.2. In the Euchdean case, Carleman estimates with hnear 
weights can be obtained from Carleman estimates with pseudocon- 
vex Carleman weights by scaling. Indeed, suppose that the following 
Carleman estimate 

|Ur(xi+x?/2£)^|| < C;^||e"(^^+^'?/2=)AM|| 2n 

II IIl7^(R") ~ II "L7JT3{Rn)' 

holds for all e < Eq and all u G C^{K), then applying this estimate to 
Ufj_ = u{ix ■) with /i > 1 and u G C^i^K)^ one gets 



iil7i=^(R") ~ II iil^^+3(R") 



Choosing yU = ^/r, and using the fact that e^^i/^^ ~ on K, one gets 
the Carleman estimate 

lle^^lMll 2n < C^' JIc'^'^IAmII 2n 

for all u G C^(iir). However, in the anisotropic case, one has to find 
another way. 

To prepare for the proof of Proposition 12.11 consider the Laplace- 
Beltrami operator on A^, 

and the corresponding conjugated operator (by the limiting Carleman 
weight Xi) 

(2.1) e^^iPe-^^i = dl^ - 2Td,, + + A,„. 



We denote by Aq = < Ai < A2 < • • • the sequence of eigenvalues of 
— Agg on Mo and (V'j)j>o the corresponding sequence of eigenf unctions 
forming an orthonormal basis of L^(Mo), 

We denote by vTj : L^(Mo) — )■ L'^{Mo),u {u,ipj)ipj the projection on 
the hnear space spanned by the eigenfunction ipj so that 

00 00 

j=0 j=0 

and by 

u{j) = / uipjdVg^ 

J Mo 

the corresponding Fourier coefficients of a function u on Mq . We define 
the spectral clusters as 

(2.2) Xk= Yl ^^N- 

k<^yXj<k+l 

Note that these are projection operators, xl — Xk, and they constitute 
a decomposition of the identity 



(2.3) ld = Yxk- 



k=0 

We end this paragraph by recalling the spectral cluster estimates of 
Sogge [211 [26j that we will need: 

(2.4) \\xku\\^^^^^^^^ < C(l + k)^-^\\u\\L2^Mo), 

\\XkU\\LHM„) < C(l + kf^-^\\u\\^^^^^^^. 

The first estimate is given in ^26^ Corollary 5.1.2] and the second one 
is a consequence of the first one by duality. 

Proof of Proposition \2.1\ Recall that our main goal is to prove 

(2.5) IImII _2n < Ct\\ I II 2n 

^ ' " "l^^(RxMo) ~ "l7^(RxA/o) 

when M G C^(/ X Mo) and 

(2.6) Dl^u + 2irD,,u - A^^u - t\ = f 

with = —idx^. The inverse operator in (12. 6p is actually easy to 
write down, as was done in [13]. The same procedure appears in [9] 



and |26l Section 5.1]. Writing / = Ylf=Q'^jf similarly for m, the 
equation formally becomes 

for Xi on the real line and for j > 0. The symbol of the operator on 
the left is + 2zr^i — + Xj, and this is always nonzero provided that 
7^ Xj for all j. Thus, if 

^ Spec(-A,J, 
an inverse operator may be obtained as 



G'r/(a;i,x') = ^ / mr{xi - yi, ^/X^) njf{yi,x') dy 



j=0 

where 

1 /""^ e^*^ 
27r + 2irr7 - t'^ + fi 



mr{t,fx) = —l r— T^dr/, ^ > 0. 



The operator is the same as Gr in [13i Section 4], except that 
in the present setting {ipj} is a basis of L^(Mo) on a compact man- 
ifold (Mo,(7o) without boundary instead of being a basis of Dirichlet 
eigenf unctions on a compact manifold with boundary. Let 

L^(R X Mo) = {fe LL(R X Mo) ; (1 + xD'^'f G L^R x Mo)} 

and let -f^|(R x Mo) by the corresponding Sobolev spaces. The proof 
of Proposition 4.1] goes through for Gr without changes and shows 
that for any fixed 6 > 1/2, if |r| > 1 and ^ Spec(— A^q) then the 
equation 

e™n-A9)e"™'^ = / 

has a unique solution v = Grf G H^g(Rx Mq) for any / G ^^(Rx Mo). 
Further, v G ^^^^(R x Mo) and 

\Mh-_^{-RxMo) < '^okl'"i/||L2(RxAfo)' < S < 2. 

We define 

Grf{Xi,x') = xiXi)Grf{Xi,x') 

with X ^ C'^(R) which equals 1 on I. It is then clear that all the 
statements in the proposition except for the estimate follow from 
the results for Gr explained above. 

It is sufficient to prove the estimate in the case where r > 4 and 

^ Spec(— Agg). We first record a lemma. 
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Lemma 2.3. // r > 0, /x > 0, r 7^ /i and t G R then 



Besides, if t > then 

|m^(t,0)| < |t|e-^l*l. 
Proof. This follows by writing 



{iri - {t + fi)){iri - {t - fi)) 2/x 
and by noting that for a > 



ir] - {t + n) ir] - {t - n) 



^-1 



1 



irj + a 



and similarly for a < 0. 
Furthermore we have 



^-1 



it) 



it) 



0, t < 
e""*, t > 0, 



te-^l*l, t<0 
0, t>0, 



□ 



and this concludes the proof of the lemma 

From the decomposition f l2.3p . the spectral cluster estimate (12 ■4p . 
and the fact that spectral clusters are projections {xl — Xk), we get 
the following string of estimates 



(2.7) 



' "L7I^(Mo) 



2n 



< ^(1 + /c)^""||XfcW|U2(Mo)- 



k=0 



Since 



\XkU\\L2(Mo) 



k<WXj<k+l 



if we apply the inequality (12. 7p to m = Grfixi, •), we get for almost 
every Xi G / 



ICr f (a^i7 Oil ^ 



<J2(l + k) 



1—1 

2 n 



fc=0 



X 



A;<,/A,<fc+1 



"^r(a;i -yi, \^j) fiyi,j)dy 
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By Minkowski's inequality, we have 



fc=0 



X 



( \'^r{xl-yl,^/\J)f{yl,j) j dyi 

k<^/>^<k+l 



and since 



mr{xi - yi, \^j)f{yi,3) 

k<^y\j<k+l 

< sup |m^(xi - yi, a/A~)|^ 5Z I -^O I 

k<^<k+l k<^<k+l 



< sup |m^(xi - v^A")! X ||xfc/(yi,-)lli2(Mo) 



using once again the spectral cluster estimate <\2Ah . we finally get 

oo 

(2.8) Oil -,^,^< 5^(1 + A:)i-i 



'L^(M„) 

fc=0 



X 



sup |m^(xi - yi, VA~)| X ||/(yi,-)||^-^ d?/i. 



Using Lemma [2.3[ we estimate 

{e"(''~^^l*l when t < k 

1 when k < T < k + 1 

^-{r-k-m ^Yien T>k + 1 

with k > 0. (Note that when k = 0, a majorant is e~^^/^-"*' for r > 4). 
This allows us to estimate the series 

oo 

^(1 + A;)"'^"" sup \mr{t,y^)\ 

k=0 k<y/x]<k+l 

2 iM.l 2 



< /t-^e-(^-'^-i)l*l+r-^+ 5^ /t-^e-^'^-^^l^l+e-^-/^)!*! 

l<fc<T— 2 /i;>r+l 

/r— 2 POO 
^-f e-("-"-2)l*l dr + 1 + y r-^e-(""")l*l dr. 
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By an obvious change of variables we have 



T— 2 noo 

r-f e"^^-'-^)!*! dr + j r-ie-^^"^)'*! dr 

•» 1 /»oo 



whence 



(2.9) + sup v^)| <l + |t| 

k = k<y/>^<k + l 

From the estimates (12. 8p and (12.91) . we obtain 

\\Grf(Xl,-)\\ _2n^ ^ 

II ^ ^iIl^i=^{Mo) 

/■oo 

< 



/OO 

/oo 



2n 



and we conclude using the Hardy-Littlewood-Sobolev inequality 

WGrfW 2" < 11/11 2" 

This completes the proof of Propositio ri2.1[ □ 



3. Complex geometrical optics 

In this section we will construct the complex geometrical optics so- 
lutions that will be used to recover an L"/^ potential. Throughout the 
section, let (M, g) be a compact manifold with smooth boundary, and 
let {M,g) m {T,g) m {f,g) where T = R x Mq, T = R x Mq, and 
g = eO) go, and (Mq, go) d (Mq, go) are two {n — l)-dimensional simple 
manifolds. We also assume that n > 3. 

First we state a consequence of Proposition 12.11 for the manifold M 
(this follows easily by embedding (Mo,5'o) in some compact manifold 
without boundary and using suitable restrictions and extensions by 
zero) . 

Proposition 3.1. For |t| > 4 outside a countable set, there is a linear 
operator Gr : L^{M) H^{M) such that 

e™^(-Ag)e-™-iG,t; = v forve L\M), 
G,e™i(-A3)e-™i?; = v forve Co"^(M"*). 
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This operator satisfies 

Co 

||G'^/||l2{M) < 1^11/11x2(4^), 
\\Grf\\H^M) < ^011/11x2(4/), 

\\GJ\\ 2„ <Co|l/i| 2n 

where Cq is independent of r. 

Let us first construct the required complex geometrical optics solu- 
tions for the case where no potential is present. This is analogous to 
[U Proposition 5.1] for g = 0. 

Proposition 3.2. Let uj G Mq\Mq be a fixed point, /et A G R be fixed, 
and let b G C°°(S'"~^) be a function. Write x = {xi,r,6) where {r,6) 
are polar normal coordinates with center u in (Mo,(7o). For |r| suffi- 
ciently large outside a countable set, there exists Uq G H^{M) satisfying 



—AgUo = in M, 

Uq = e"™i(e"*"'|^| '''e-'~-~ '0[t>) +roj 



where tq satisfies 

kllko||L2(Af) + \\ro\\m{M) + \\ro\\^^^^^^ < 1. 

Proof. The claim follows if one can find vq satisfying 

e™H-A>"™Vo = / 
with the required norm estimates, where 

/ = e™iAj,e-™i(e-^^^e^^(^i+^'^)6(^)). 
It is enough to take tq = Grf and to note that 

U.(M) = ||A,(e^"(^^+'^)&W)IU2(./)<l. 



The LF' and estimates follow from Proposition 13.11 The esti- 
mate follows from the estimate and Sobolev embedding, or alter- 

2n 2n 

natively from the L"+2 — )■ L"-2 estimate for Gr- □ 

We next consider the case with a potential q G L"'/^(M), and try to 
find a solution to (— + q)u = in M of the form 

u = uo + e"™Vi. 

Since — Ag-ug = 0, the function ri should satisfy 

(3.1) e™^ (-Ag + g)e-"^Vi = -ge^^Mo- 
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2n_ 
2 



Since q is only in L"/^(M), here we need to use the L"+^ — )■ L" 
estimates for Gr- We follow the argument of Lavine and Nachman 
[T7] . It will be convenient to symmetrize the situation slightly. Later 
on, the L" functions rrij in the next lemma are chosen to be essentially 

Lemma 3.3. Let mi,m2 G L"^{M) be two fixed functions. Then for 
|t| > To outside a countable set, 

(3.2) \\miGrm2f\\L'^ < Co||mi||Ln||m2||Ln||/||i2. 
Further, 

(3.3) \\miGrm2f\\L2{M)^L2{M) 

as \t\ — )■ oo. 

Proof. The Holder inequality and Proposition 13.11 imply that 

||miG^m2/||L2 < ||mi||in||G^m2/|l ^ < Co||mi||in||m2/|l ^ 



< C'o||"^i||l"||"^2||l"||/||l2. 

Let £ > and decompose nij = nij + m^- where m^- G L°°(M), 

II'^j'IIl" < ll'^jlU" 

and \\m)\\L^ < — -—^ . . — -. 

3Comax(||mi||Ln, ||m2||L"j 

(One can take for instance m^- = f^jXHnijiKfi} for large enough fi.) It 
follows from the estimates for Gr and (13.21) that 

\\miGrm2f\\L2 < \\m\G-rmlf\\L2 + \\m\Gr'my\\L2 + \\m\Grm2f\\L'^ 

( Co\\m\\\L^\\m\\\L^ e e\ 

- \ M s + sj I'^l'^^- 

The last expression is bounded by £:||/||£,2 if |r| is sufficiently large. 
This proves (13. 3p . □ 



We now finish the construction of complex geometrical optics solu- 
tions. 

Proposition 3.4. Assume that q G L"/2(M). Let uj e Mq\Mo be a 
fixed point, /ei A G R be fixed, and let h G C°°{S^~'^) be a function. 
Write X = {xi, r, 6) where (r, 6) are polar normal coordinates with cen- 
ter u in (Mo,(7o). For \t\ sufficiently large outside a countable set, there 
exists a solution u G H^{M) of (— + q)u = in M of the form 

u = e-™i(e"'^1c/r^^V^^'''"^'''^&(^) + f) 
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where f satisfies 

\\f\\ 2n < 1, ll^llLsrA-n — !■ as Irl — )■ oo. 

Proof. As explained above, we let uq be the harmonic function given in 
Proposition l3.2[ and look for a solution of the form u = Mo+e~™^ri. We 
write q{x) = \q{x)\e^°'^^'' = \q{x)\^^'^m(x) where m(x) = \q(x)\^^'^e^°'^^\ 
Then |g|^/^,m G L"{M) with L" norms equal to Hf^H^i^. 

We obtain a solution u provided that (13. ip holds. Trying ri in the 
form ri = G^-lgl-'^/^f , we see that v should satisfy 

(Id + mGr\q\^^'^)v = -me'''^UQ. 

By Lemma 13.31 for |r| sufficiently large one has HmGrlgl^^^HLZ^LZ < 
1/2. One then obtains a solution 



-(Id + mGrW^)~^{me^''^uo) 

1. Consequently 



Since ||me™^Mo||L2 < H'riULn ||e™^Mo|| 2n < 1, it follows that ||f ||l2 < 



llrill 2n < Cnll lO'l^^^t'll 2n <1. 

II ^iIltt^ — ^iii^i iIlht2 ~ 

Now u is of the form given in the statement of the proposition, provided 
that 

f = ro + ri . 

This remainder term satisfies ||f|| 2n < 1. 

II iilh=^ ~ 

To study 11^11^2 we fix e > and make a decomposition = s'^ + s^ 
where s« e L°^{M), Ws^l-^ < Ugll^^a, and \\s^\\Ln < e. Then 



kl||L2 < ||G'T-S''t>||j;^2 + Cl ||G'T-s''f II 2n 



, / Co||s^||l°° , ^ ^ II bii \ II II 



Choosing |r| sufficiently large, we see that ||ri||i2 < e for |r| large. 
Since also ||^'o||l2(m) ^ I'^l^^j it follows that ||r||2,2 — )■ as |r| — t- 00. 

Finally, to prove that u e H^{M), it is enough to consider a compact 
manifold (M, g) which is slightly larger than (M, g) and extend q by 
zero outside M, and to perform the above construction of solutions in 
M. One obtains a solution u G L'^{M) C L'^[M), and A^m = qu E 
L"+2 (M) C H~^{M) by Sobolev embedding. Elliptic regularity implies 
that M e i//„,(M'°*), thus also u G H\M). □ 
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4. Recovering the potential 
We are now ready to give the proof of the main uniqueness result. 

Proof of Theorem \l.l[ Assume, as before, that {M,g) d (T,g) d (T,g) 

where T = R x Mq, T = R x Mq, and {Mo,go) d {Mo, go) are two 
(n — l)-dimensional simple manifolds where n > 3. Also assume that 
g = e®go. 

From the assumption A^^^ = A^^^, writing q = qi — q2, we know 
from Lemma [A. II that 

(4.1) / qu,U2dVg = 

J M 

where Mi,U2 € H^{M) are solutions of (— + qi)ui = in M and 
(— + q2)u2 = in M. By Proposition 13. 4[ for r sufficiently large 
outside a countable set there exist solutions Uj of the form 

Ml = e-"("^+''^)(|^r^/V^("^+''")6(^) +ri), 

U2 = e^^^^^^'^'\\g\-'/' + r2). 

Here A is a fixed real number, b G C°^(S'"^^) is a fixed function, and 
X = {xi,r,6) are coordinates in T where {r,6) are polar normal coor- 
dinates in (Mo, go) with center at a fixed point u G Mq \ Mq. Also, the 
remainder terms satisfy 

as r — )■ oo. 

Inserting the solutions in fl4.ip and noting that dV^ = |(7|^/^dxi dr d9, 
we obtain that 

(4.2) [ qe'^^'''^''^b{e)dxidrde = I g(air2 + aan + nrs) d1/ 

where ai,a2 are smooth functions in M independent of r. We show 
that the right hand side converges to as r — )■ oo. To do this, fix 
£ > and write q = q^ + (^ where q^ G L°°(M), ||q'^||Ln/2 < ||g'||Ln/2, and 
||?''||l"/2 < s. The right hand side of (14.21) is bounded by 




g(air2 + a2ri + rir2) dV 



^ Ik IU°°(lki||L2 + ||r2|U2 + ||ri||i2||r2||L2) 

+ lkllL"/2(lkl|L^ + 1^2 IL^ + ||ri|| 2n llrall 2n ). 
irJ II-IJ "L"-2 "L"-2 "L"-2 "L"-2 

Using the bounds for r^, if r is sufficiently large then the last quantity 
is < e. This shows that the right hand side of (14. 2p goes to as r — t- oo. 
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Extend q by zero into T and interpret the left hand side of (14 .2^ as 
an integral over T. Taking the limit as r — t- oo, we obtain that 

/OO /'OO p 
/ / g(xi, r, ey^^'''+'%{e) dxi dr d^ = 0. 

This statement is true for all choices of G Mo\Mo, for all real numbers 
A, and for all functions b G C~(S'"-2). Since q G ^^(M), Fubini's 
theorem shows that q{- ,r, 9) is in ^^^(R) for a.e. (r, 9). Consequently 

r poo 

(4.3) / / /A(r,^)e-^''6(^)drde = 

where fx G L}{Mq) is the function given by 



/oo 
e*^^ig(xi,r,^) dxi. 
'OO 



If |A| is sufficiently small, it follows from Lemma l5.ll below that the 
vanishing of the integrals (14. 3 p for all choices u and h implies that 
fx = 0. Since q{- ,r, 9) is a compactly supported function in i^^(R) for 
a.e. (r, 9), the Paley- Wiener theorem shows that q{- ,r,9) = for such 
{r,9). Consequently qi = q2- □ 

5. Attenuated ray transform 

It remains to show the following lemma which was used in the proof 
of Theorem 11.11 

Lemma 5.1. Let (Mq^qq) be an {n — 1)- dimensional simple manifold, 
and let f G L}{Mq). Consider the integrals 

f{r,9)e-^''b{9) dr d^ 







where {r,9) are polar normal coordinates in {Mq^qq) centered at some 
(jj G OMq, and t{u, 9) is the time when the geodesic r i-> (r, 9) exits Mq. 
// |A| is sufficiently small, and if these integrals vanish for all u G dM^ 
and all b G C^{S''-'^), then / = 0. 

The last result is related to the vanishing of the attenuated geodesic 
ray transform of the function / on Mq. For the following facts see [3], 
[20j . [22j . To define the ray transform, we consider the unit sphere 
bundle 

SMo = \J S,, ^, = {(a;,OGT,Mo; 1^1 = 1}. 
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This manifold has boundary d{SMo) = {(x, ^) G SMq; x G dMo} 
which is the union of the sets of inward and outward pointing vectors, 

d±{SMo) = {{x,0 e SMo;±{^,u) < O}. 

Here u is the outer unit normal vector to OMq. Note that 9+(S'Mo) is 
a manifold whose boundary consists of all the tangential directions 
e d{SMo); = 0}. Thus the space C^{{d+{SMo)y''') 

contains all smooth functions on 9+(5'Mo) vanishing near tangential 
directions. 

Denote by t h-> 7(^,0;,,^) the unit speed geodesic starting at x in 
direction ^, and let t{x,^) be the time when this geodesic exits Mq. 
Since {MQ,go) is simple, t{x,^) is finite for each (x,^) G SMq. We also 
write {pt{x,^) = {■y{t, x, ^) , 'j{t, x, ^)) for the geodesic flow. 

The geodesic ray transform, with constant attenuation —A, acts on 
functions on Mq by 

Txf{x,0= r^'^^ f{lit,x,0)e-^'dt, {x,Oed+iSMo). 
Jo 

In Lemma [5.11 if / were a continuous function, one could choose b{6) 
to approximate a delta function at fixed angles 6 and obtain that 

T{uj,e) 

f{r, e)e-^' dr = 

for any u G OMq and any 6 G S*""^. Since (r, 6) are polar normal 
coordinates the curves r {j-,6) are geodesies in {Mq^qq), and this 
would imply that 

Ta/(x,0 = for all G 9+(^Mo). 

One has the following injectivity result from [4, Theorem 7.1]. (If 
Mq is two-dimensional the smallness assumption on the attenuation 
coefficient was recently removed in [2T].) 

Proposition 5.2. Let [MQ^go) he a simple manifold. There exists 
£ > such that if X is a real number with |A| < e and if f & C°°{M), 
then the condition Txf{x,^) = for all {x,^) G 9+(S'Mo) implies that 

f = o. 

The previous argument together with Proposition 15.21 proves Lemma 
15.11 for smooth /. However, this requires well defined restrictions of 
/ to all geodesies and it is not obvious how to do this when / G L^. 
We circumvent this problem by using duality and the ellipticity of the 
normal operator T^Tx- 
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We will need a few facts. Below we write 
for h G C°°((9+(5Mo)), and 



ih,h)Li{d+(SMo)) = hhij,d{d{SMo)) 

Jd+{SMo) 

where fi{x,C,) = — z/(a;)) and dA^ is the natural Riemannian volume 
form on a manifold N. 

Lemma 5.3. (Santald formula) If F : SMq — > R is continuous then 
Fd{SMo) 

SMo 



ld+(SMo) Jo 

Proof. See P| Lemma A. 8]. □ 
Lemma 5.4. If f E C°^(Mo) and h G Co^((9+(5Mo))"*) then 

{Txf, /i)L2(a+(5A/o)) = {f,T^h)L2i^Mo) 

where 

nh{x)= I e-^^(^'-«)/i^(x,0d5',.(0, xgMo. 

J Sx 

Proof. By the Santalo formula 

(Ta/, h)Ll{d+iSMo)) 



d+{SMo) Jo 



e-^V(7(t,a:,0)Vdtd(a(5Mo)) 
e-^Vlv't 0)hi,Mx, dt d(5(^Mo)) 

a+(5Afo) Jo 

e-^^(^'-«/(x)/i^(x,Od(^Mo) 

/(x)f / e-^-(^'-«)/i4x,0d5.(0) dV{x). 

'Mo \JS^ / 

This proves Lemma [5 .41 □ 

Lemma 5.5. T^Tx is a self-adjoint elliptic pseudodifferential operator 
of order —1 in Mq"* . 

Proof. This is contained in [6l Proposition 2], but for completeness we 
also include a proof in Appendix [Bl □ 



19 



Proof of Lemma \5. il The first step is to extend (Mo,(7o) to a slightly 
larger simple manifold and to extend / by zero. In this way we can 
assume that / is compactly supported in M™*. 

We let h also depend on u and change notations to write the assump- 
tion in the lemma in the form 

'^^'^^ e-^V(7(t, X, 0)b{x, dt dS,{0 = 

's^ Jo 

for all X G dMo and b G C^{{d+{SMo))"'^). Next we make the choice 
= h{x,^)fi{x,C) for h G C^Hd+iSMo))'''^) and integrate the 
last identity over dMo to obtain 

f r^"'^^ e-^V(7(t, X, 0)Hx, 0/i dt didiSMo)) = 0. 

Jd+{SMo) Jo 

We are now in the same situation as in the proof of Lemma 15.41 and 
invoking the Santalo formula implies 

f{x)T^h{x)dV{x) = 

Mo 

for all h G C^((a+(^Mo))^^*). Note that the last integral is absolutely 
convergent because / G L^(Mo), and also the previous steps are justi- 
fied by Fubini's theorem. 

It remains to choose h = T\ip for G C^(M™*) to obtain that 



f{x)TlT^^{x) dV{x) = Q. 

I Mo 

Since T^Tx is self-adjoint, we have 

{TlTxf{x))^{x)dV{x) = Q. 

Mo 

This is valid for all test functions so T^T\f = 0. By ellipticity, since 
/ was compactly supported in M™*, it follows that / G C^{Mq^^). 
One can now use the argument for smooth / given above, together 
with the injectivity result (Proposition 15. 2p . to conclude the proof that 
/ = 0. □ 

Appendix A. Wellposedness 

Here we recall the standard arguments that show wellposedness of 
the Dirichlet problem for —Ag + q on a compact oriented manifold 
{M,g) with smooth boundary and with q G L'^/^(M), n>3. Consider 
first the inhomogeneous problem for the Schrodinger equation, 

(A.l) + q)u = F in M, uIqm = 0. 
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The bilinear form related to this problem is 

B{u, v) = / {{du, dv) + quv) dV, u,v e H] (M) , 

J M 

where ( ■ , ■ ) is the complex-linear inner product of 1-forms and dV" 
is the volume form on {M,g). By the Sobolev embedding H^{M) C 

L"-^{M) and by Holder's inequality, 5 is a bounded bilinear form on 
Hl{M). Writing g = g« + where e L°°(M) and \\q^\\L^/2(^M) is 
small, we obtain from Poincare's inequality that 

B{u,u) > c\\u\\]j,(^M) - C\\u\\l2(My u e Hl{M). 

This shows that -B + C is coercive, and by the Lax-Milgram lemma, 
compact Sobolev embedding and the Fredholm theorem, the equation 
( ]A.1|) has a unique solution u G Hl{M) for any F G H^^{M) if one is 
outside a countable set of eigenvalues. 

We can now consider the Dirichlet problem 

(A.2) {-Ag + q)u = in M, uIqm = /• 

We assume that is not a Dirichlet eigenvalue, and it follows from the 
above discussion that for any / G H^/'^{dM) there is a unique solution 
u G H^{M). The DN map is formally defined as the map 



\dAI- 



f ^ d^u\ 

More precisely, if / G H^^'^{dM) we define Ag^gf weakly as the function 
in H~^/'^{dM) which satisfies 



dM 



{KgJ)hdS= I {{du,dv) +quv)dV 

J M 



where u is the unique solution of flA.2p . and v is any extension in H^{M) 



of h {v\9M = h). Then Ag^g is a bounded map H^/\dM) H-^/\dM) 
again by Holder and Sobolev embedding. 
The DN map satisfies in the weak sense 



[ {AgJ)hdS= [ fAg,ghdS. 

JdM JdM 



To see this, let u,v E H^{M) solve (— A^, + q)u = 0, u\sm = f and 
(-Ag + q)v = 0, v\gM = h. Then 

{Ag^gf)hdS= [ {{du,dv) + quv)dV 

dM ' Jm 



[ {{dV,du) +qVu)dV = [ {Ag^gh)fdS. 

Jm JdM 
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As a consequence, we have the basic integral identity used in the 
uniqueness proof. 

Lemma A.l. Ifqi,q2 e L"'/^(M) and kg = Ag^q,^, then 

/ {Qi - q2)uiU2 dV = Q 
J M 

for any Uj G H'^{M) with {-Ag + qi)ui = in M , {-Ag + 52)^2 = 
%n M. 

Proof. Follows from the computation 

= / (Ag^g^ - Ag^g2){Ul\9M)U2dS 

JdM 



{Ag,q^{Ui\oM)u2 - UxAg^g^{u2\dM)) 

dM 

and the definition of the DN maps. □ 



Appendix B. Normal operator 

The setting is the compact simple Riemannian manifold (Mo,(7o) of 
dimension n — 1. Let Tx be the attenuated ray transform as in Section 
[5l We will prove Lemma 15.51 Write 

We compute the normal operator T^Txf for / G C^(Mg'^*) 

nTxfix) 



e 



(•r{tp{x,(,)) 



JO 
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and using changes of variables we get for tlie last integral expression 

Js^ J-t{x,-0 



/ / + / 



[e 



Changing variables y — exp^{s^) shows that 



K,{x,y)fiy) dViy) 



Mo 



where 

K\{x,y) 
with 



(^Q-X<p+{x,y) ^ Q-X<P-{x,y)^ /detyo(^ 



dgo \x^y) 



detg'o(y) 



det(exp^ )'{x,y)\. 



(fi± = 2t{x, ^ gradj^ dg^{x, y)) ± dg^{x, y). 

The functions (p± are smooth away from the diagonal x = y, and their 
A;-th order derivatives behave y) ^. Note that det(exp 

stands for the Jacobian determinant of 

exp-i : Mo X Mo ^ R'^"^ 

{x,y) ^ exp-^(y). 

The kernel of the normal operator is symmetric 

Kx(y,x) = Kx{x,y) 

and the singular support of this kernel is the diagonal in Mq x Mq. 

We will now prove that the operator Px with kernel Kx is actually 
a pseudodifferential operator. The first observation in that direction is 
that in coordinates 

(B.l) dl^ix,y) = ajkix,y){x^ - yj){x'' - y^) 

with a^^{x^x) = gi^{x). Indeed the square of the distance vanishes at 
second order and its Hessian aX. x = y is twice the metric. This can be 
seen from the well known formula 

VV(y)(^,^) = ^<^(exp,i^)[^^ 
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and the fact that if \6\g^ = 1 then d'^^{expyt9,y) = t^. To prove that 
Pa is a pseudodifferential operator in \[^~^(M™*) we need to show that 
for any couple of cutoff functions (V'jX) supported in charts of M™*, 
the operator with kernel 



Kx{x,y) = ip{x)Kxix, y) a/ det go{y)x{y) 

expressed in coordinate^], is a pseudodifferential operator on R"~^ with 
symbol in S^^. Because of its form and of (IB.ip . the kernel satisfies 

(B.2) \d:d^kx{x,x~y)\<C, 



\y 



and has compact support in R"~^ x R"~^. 

Such operators are pseudodifferential operators and this can easily be 
seen in the following way: the symbol associated with such an operator 
is 

px{x,^) = j kx{x,x - y)e~^y'^ dy 

For cutoff functions ip and x whose supports don't intersect, the pre- 
vious symbol is a Schwartz function because the kernel is a smooth 
compactly supported function. So we are only interested in those sym- 
bols corresponding to kernels Kx{x,x — y) which are supported close 
to R"~^ X {0}- In that case, we use a dyadic partition of unity 



1 = E ^(2" 



with X a function supported in an annulus, to decompose the symbol 
as a sum of terms of the form 

2M(n-i) f ^^'y<^{y)Kx{x,x~2^^y)dy. 



Note that because of the compact support of the kernel, these terms 
vanish when /i is large, so we are mainly concerned with the terms 
where is less than some positive integer, say A^. Because of the be- 
haviour (lB.2p . the rescaled kernel Kx{x, x — 2^y) is uniformly bounded 
by 2^'*'^""^) as well as all its derivatives. Applying the non-stationary 
phase when |^| > 1 and 2^^ is large we get 

Ipa(x,OI< E 2^(21^1)-"^+ E 2^ 

<(l + lel)-^ 



^By a slight abuse of notations, to lighten the exposition, we don't write the 
pullback by the coordinates and think of x and y as variables in R"~^. 
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Repeating this argument for the derivatives of this function, we get 
that px is a classical symbol of order —1. 
Let us concentrate on po- we have 

So y) 

from the previous computation, we see that taking x = y in the non- 
singular factors, yields error terms whose kernel are less singular by an 
order of \x — y\, i.e. errors with symbols of order (1 + Therefore 
in terms of the principal symbol, it suffices to compute 

iIj{x)x(x) X y/detgo{x) / —— ^dy 

J \Q()ix)y-y\ 2 



Cn\9o\x)C-C\ ^iIj{x)x{x) 



Finally, these observations show that Pq has as principal symbol a 
multiple of 



\f^ix)i^k 



and since the principal symbol of P\ depends smoothly on A, it doesn't 
vanish for A small enough. This means that for A small enough, Px is 
an elliptic self-adjoint pseudodifferential operator of order —1. 
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